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ABSTRACT 


Charts of potential vorticity (g) within four 
different density layers, calculated from an eddy resolving 
global ocean general circulation model, are presented for 
the Atlantic, Pacific, and Indian Ocean Basins. Potential 
vorticity is evaluated as g = (f/p) * (Acr/Az), a formulation 
that is valid for large, slow scales and allows it to be 
readily calculated from model output. Here, f is Coriolis, 
p in situ density, <5 the potential density referred to the 
sea surface, and z the vertical coordinate. 

The model is shown to represent oceanic density 
structure with reasonable accuracy. The most significant 
departure from observations are the steep density gradients 
on the northern boundary of the Antarctic Circumpolar 
Current (ACC). This may be due to model convective 
adjustment parameterizations and the extent of the smoothing 
in the cataloged data, which may be greater in the southern 
oceans due to relatively scarce observations. Model derived 
values of potential vorticity correspond well with those 
calculated from observations. In the middle and lower 
layers, the three basins tend toward homogenization, 
representing unventilated layers. "Gyres" of g in the Nil: 


with their northern boundaries anchored in the outcrop 
region, illustrate the ventilation of these layers. 
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I. INTRODUCTION 


Several theories attempt to predict the horizontal and 
vertical structure of the thermocline (e.g., Rhines and 
Young, 1982; Luyten et al., 1983; Pedlosky and Young, 

1983; Pedlosky, 1983). These theories all revolve around 
the role of potential vorticity (g) . Potential vorticity 
is a calculated oceanic property that is conservative in the 
absence of processes which change water parcel 
characteristics, such as turbulent mixing or surface 
forcing. It can be considered a measure similar to potential 
angular momentum induced by Coriolis and other forces 
affecting water properties, as reflected primarily in 
changes in density gradients. 

The Ertel potential vorticity is the total vorticity 
(the sum of Coriolis and relative vorticity) times the 
gradient of density, divided by the density. For most of 
the ocean, the strongest density gradients are in the 
vertical, so potential vorticity is approximately the total 
vorticity times the vertical gradient of density, divided by 
density. 

A final approximation is that the relative vorticity is 
small compared to Coriolis (assuming a small Rossby number) 
so relative vorticity can be eliminated. Calculations of 
relative vorticity from model output support this 
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approximation. McDowell, et al. (1982) defined potential 
vorticity as g = (f/p)*d(p)/dz. Where f is Coriolis, p in 
situ density, and z is the vertical coordinate. In 
isopycnal coordinates, the density gradient is represented 
by the distance between density surfaces rather than the 
density change between pressure (depth) surfaces. For large 
scale, quasi-geostrophic flows, relative vorticity is much 
smaller than the vortex stretching part of g, and so, to 
predict the field of potential vorticity is to predict the 
vertical separation between density surfaces. Determining 
the vertical structure of the global thermocline hinges 
around this problem. 

The objective of this study is to calculate model 
fields of potential vorticity from the Semtner and Chervin 
(1988, 1992) global ocean circulation model and compare them 
to observed potential vorticity fields. This will indicate 
the character of the model ocean's thermocline and the 
accuracy of its portrayal of the observed thermocline, as 
determined from the Levitus (1982) data. The intention is 
to see where the model is successful in correctly predicting 
potential vorticity distribution and to see where it fails. 
In particular, model results will be compared with the work 
of Keffer (1985) and his treatment of the "directly 
ventilated" theory of Luyten et al . (1983) and the 
"unventilated" theory of Rhines and Young (1982a). While 
these theories predict fluid properties as they propagate 
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along isopycnal layers, they differ in the degree to which 
they consider the layer to be isolated from surface boundary 
conditions. 

The degree of exposure to surface boundary conditions 
and forcing will determine the characteristics of a layer 
within the subtropical gyre thermocline. The potential 
vorticity of fluid particles recirculating within ventilated 
layers is reset with each circuit around the gyre, as the 
particles pass through the outcrop zone. This outcropping 
of density layers allows deep vertical motion to be forced 
by Ekman pumping, resulting in deep ventilation, which 
advects to latitudes far south of the outcrop where the 
layer is shielded from direct wind influence. Particles 
within unventilated layers remain isolated from surface 
conditions long enough to approach a common homogenized 
state, sharing their properties with neighboring particles. 

The two theories each have features they will predict 
and are not in direct competition with each other. 

Generally, the theories complement each other. In addition 
to the comparison with Keffer's (1985) results, model fields 
will be compared with an earlier study of North Atlantic 
potential vorticity performed by McDowell, et al. (1982). 

Section II is a brief review of potential vorticity 
theories, while in Section III, observations of potential 
vorticity are discussed. In Section IV the data base and 
potential vorticity algorithm are discussed. The actual 
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model output is introduced and described in Section V and 
model results are compared with observed potential vorticity 
fields. Conclusions follow in Section VI. 
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II. THEORIES OF POTENTIAL VORTICITY 


Theories depend on whether layers within the 
thermocline are considered "ventilated," (having an outcrop 
at the surface) or “unventilated" (no surface outcrop). 

Note that a layer which outcrops may be poorly ventilated. 
This can happen if the Ekman pumping at the outcrop is weak 
or if additional dynamics act to isolate the layer from the 
boundary conditions at the outcrop. This latter point is 
the assumption that serves to separate the theory of Luyten 
et al. (1983) from the theory of Rhines and Young (1982a). 

A layer is considered unventilated if the fluid recirculates 
several times within the layer before having its properties 
reset at the surface outcrop. Conversely, a well ventilated 
layer has its properties reset with essentially every 
circuit around the gyre. This differs from the geochemical 
definition of ventilation. To avoid implying something 
about a layer's oxygen saturation, layers are referred to as 
either directly ventilated or dynamically isolated from 
surface boundary conditions. The latter may also be 
referred to as unventilated, indicating a layer that does 
not outcrop. 

Next, the theories of Rhines and Young (1982a) and 
Luyten et al. (1983) are briefly reviewed. For greater 
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detail, the reader is referred to the original manuscripts. 
In a weakly driven layer that does not outcrop, the 
potential vorticity contours are dominated by f = f c + fiy 
and the lines of constant g will run east-west, 
intersecting meridional boundaries. Only weak flow is 
possible if the fluid is to preserve its density and 
potential vorticity (see Rooth et al., 1978). Any 
meridional transport is confined to the uppermost layer 
where the forcing is strong enough to change the fluid's 
potential vorticity and drive fluid particles across lines 
of constant g. To allow recirculation, the particles must 
first forget this new value of g and return to the old. 

This is the role of the western boundary current of 
traditional linear models where dissipation removes the 
energy added by the forcing. The net result is that in a 
weakly driven, linear ocean, circulation is only allowed 
where forcing or dissipation is strong: in the upper layer. 
In deeper layers, the g contours are blocked by meridional 
coastlines and the forcing is not there to drive the flow 
across them. 

The vertical structure of velocity is more complicated, 
implying ocean circulation is either nonlinear and strongly 
driven or, perhaps, diffusive. Modification of the g 
contours by the flow itself can allow the interior fluid to 
move meridionally. In the theory of Rhines and Young 
(1982), the contours of g close on themselves if the 
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forcing is strong enough, forming rings coinciding with the 
streamlines. As a result, fluid circulates freely around 
the gyre while preserving its potential vorticity. This 
allows the western boundary current to be inertial vice 
dissipative. 

In a following paper, Rhines and Young (1982b) applied 
an extension of the Prandtl-Batchelor theorem suggesting g 
should be well-mixed within these closed rings. For this to 
occur, there must be no sources or sinks of q within the 
closed contours, the density layer must be dynamically 
isolated or have a large recirculation index. The 
recirculation index is a measure of how many times a 
particle circulates around the gyre before having its 
potential vorticity or density reset. Therefore, g 
homogenizacion is only likely to happen where the density 
surfaces are deep enough to be isolated from surface 
forcing, but not so deep as to avoid the influence of the 
indirect eddy forcing (the u'g' terms) allowing the fluid 
to spin up and close contours. 

An alternative theory of the structure of the 
thermocline, valid for small recirculation indices, was 
developed by Luyten et al . (1983). In this theory, the 
thermocline is regarded as ventilated. That is, fluid 
properties are pumped down from an outcrop along (g, a e ) 
contours via Ekman pumping, deep into the interior of the 
subtropical gyre, while preserving density and vorticity. 
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The value of q on a density surface is then its surface 
value. Knowing q at the surface, or the surface density 
and an Ekman pumping rate, will completely specify the 
interior values of q throughout the subtropical gyre, where 
the q contours can be traced back to the surface. In such 
a ventilated flow, fluid particles do not tend to 
recirculate before having their potential vorticity reset 
(they have a small recirculation index) but are instead 
pumped directly into the thermocline. 
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III. OBSERVATIONS OF POTENTIAL VORTICITY 


Over the past two decades, several papers have used 
potential vorticity to discuss circulation patterns and flow 
paths. Behringer (1972) was likely the first to draw charts 
of g on isopycnal surfaces. Two charts of the South 
Atlantic were drawn (Figures 1 and 2). The shallower 
surface, <T e = 30.75 (referenced to 750 db), showed a broad 
homogenized region of q in the center of the subtropical 
gyre. The deeper surface, <T e = 32.1 (1000 db reference), 
also showed homogenization, but the charts hint at the 
intrusion of high g waters from the Indian Ocean. 

Behringer (1972) used a limited data set. With the more 
data rich Levitus (1982) climatological data used by Keffer 
(1985), these features survive but take on a new light, as 
will be discussed in Section V. 

Coats (1981) determined absolute flow vectors in the 
North Pacific (35° N, 155° W) using a variation of the "f}- 
spiral" method, minimizing the geostrophic noise. He found 
that in the 400-1300 m range, there were no gradients of 
potential vorticity on density surfaces in this area. Coats 
(1983) found large regions of uniform potential vorticity in 
the South Pacific in the shallow subtropical gyre. 

McCartney (1982) used the conservative property of g 
to trace subpolar "mode" deep water into the subtropical 
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gyres of the world ocean. Mode water is characterized by a 
“pycnostad": a thick density layer usually formed by deep 
wintertime convection. This layer represents a source of 
low potential vorticity, which can be traced along density 
surfaces. Two such sources of low g influence the North 
Atlantic subtropical gyre. One is the well-known 18° water 
(Worthington, 1959), formed in the intense recirculation 
regions of the Gulf Stream. The other is formed in the 
subpolar gyre and enters the subtropical gyre by passing 
under the North Atlantic Current on deeper density surfaces 
(McCartney, 1982; McCartney and Talley, 1982). Both are 
described in further detail in Section V. 

McDowell et cl. (1982) drew maps of the potential 
vorticity field of the North Atlantic in several different 
density layers and interpreted them as flow paths. As 
density increased, the transition from ventilated to 
dynamically isolated was characterized by vast regions of 
nearly uniform potential vorticity (Figure 3), in agreement 
with the theory of Rhines and Young (1982b). At still 
deeper layers, the circulation had weakened enough that 
ventilation could again take hold and the layers showed g 
contours reaching up into the high latitude outcrop windows, 
allowing easy north-south movement of fluid particles 
(Figure 4). 

Keffer (1985) analyzed the digitized climatological 
atlas of Levitus (1982) . His work shows the North Atlantic 
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thermocline to be highly anomalous. The extensive 
outcropping of its density surfaces equatorward of the zero 
wind-stress curl line, and the resulting strong downward 
Ekman pumping, allows the surface boundary conditions to 
control the interior thermocline structure as suggested by 
Luyten et al. (1982). The other oceans show a greater 
tendency toward homogenization, even at great depth, 
suggesting relatively complete dynamical isolation. Only in 
the shallowest layers, where fluid particles are exposed to 
direct atmospheric forcing is there a tendency for outcrop 
properties to propagate into the interior. In these 
regions, winter-time stratification properties are advected 
into the interior as tongues of high or low g. These 
results are described in further detail in Section V. 
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IV. DATA BASE AND METHODS 


The data base is model output from the Semtner and 
Chervin eddy resolving global circulation model (Semtner and 
Chervin 1992) averaged over the last five years of the model 
run. An earlier 20-year spin up (Semtner and Chervin 1988) 
was extended for 12.5 additional years: the first 2.5 years 
with continued annual mean forcing and the final 10.0 years 
with climatological monthly forcing using European Center 
for Medium-Range Weather Forecasting (ECMWF) winds 
(Trenberth et al, 1989) . The 0.5°, 20-level grid model 
produced remarkably realistic three dimensional flow fields, 
although high latitude eddy resolution was marginal. 

Seasonal forcing improved the simulation, most notably in 
the tropics and high northern latitudes. Mid-latitude gyre 
circulations, western boundary currents, zonal equatorial 
flows and the Antarctic Circumpolar Current (ACC) exhibit 
mean and eddy characteristics similar to those observed. 

Mean flow eddy intensification of the ACC and separated 
boundary jets are also present. 

The model data base was analyzed using the National 
Oceanic and Atmospheric Administration (NOAA)/Pacific Marine 
Engineering Laboratory {PMEL) software tool FERRET. FERRET 
is an interactive program designed for use with very large, 
well-ordered data sets, allowing the user to work with 
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explicitly stored variables, abstract mathematical 
functions, mathematical transformations, and combinations of 
disk variables. Potential density (O) was calculated using 
a FERRET journal file adapted from the simplified equation 
of state used in the model and the hydrostatic balance. 
Potential vorticity was calculated as q=f*d(a)/(dz*p)* 10 13 , 
where dz is the vertical distance between two density 
surfaces differing in density by d(O). 

The entire procedure was followed for the four density 
layers listed in Table 1. These density layers were chosen 
to match those of Keffer (1985) representing depths within 
the directly forced (layer A), homogenized (B), and isolated 
(C and D) regions. 


TABLE 1. DENSITY LAYERS USED IN THE POTENTIAL VORTICITY 
CHARTS. 





Midlayer 

density 

Density 

Layer 

Ocean 

Figure 

(O e ) 

Interval (a) 

A 

Atlantic 

5 

26.40 

26.30-26.50 


Pacific 

6 

26.15 

26.05-26.25 


Indian 

7 

26.40 

26.30-26.50 

B 

Atlantic 

14 

26.75 

26.50-27.00 


Pacific 

15 

26.50 

26.25-26.75 


Indian 

16 

26.75 

26.50-27.00 

C 

Atlantic 

22 

27.15 

27.00-27.30 


Pacific 

23 




Indian 

24 



D 

Atlantic 

31 

27.40 

27.30-27.50 
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V. RESULTS 


A. UPPER THERMOCLINE (LAYER A) 

The pressure of the central density surface (o 6 = 26.4 
in the Atlantic and Indian Oceans, 26.15 in the Pacific 
Ocean) of the upper-thermocline layer is shown in Figures 5- 
7. Of note are the steep slopes at approximately 40° S 
across the South Atlantic and Indian Oceans. These near¬ 
vertical slopes are primarily due to convective adjustment 
in the model. There is some indication from observations 
that convective overturning produces similar gradients in 
these areas. However, observations remain too sparse to 
reach definitive conclusions. Also of interest are the 
shallow equatorward shelves, or plateaus, of density in the 
eastern North Atlantic and southern Indian Oceans. These 
are indicative of tongues of cold water extending toward the 
equator. In the Atlantic and Pacific Oceans, the model 
density surfaces follow the observed shape: surfacing at the 
outcrop regions, reaching their greatest depths near 30° N 
and S, then shoaling toward the surface at the equator. 

Figures 8-10 illustrate the corresponding charts of 
potential vorticity. Water parcels in this layer pass 
through the outcrop zone with nearly every circuit, having 
their properties strongly influenced before they are 
subducted back into the interior. This "ventilation" of the 
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thermocline is well illustrated in the Pacific Ocean by 
loops of relatively high potential vorticity anchored in the 
outcrop region (Figure 9). In the Atlantic Ocean (Figure 
8), the presence of Worthington’s (1959) convectively 
produced "18° water" is not as well defined in the model 
output as in Keffer's (1985) analysis of Levitus* (1982) 
data (Figure 11). 

The relative maximum in potential vorticity in the 
vicinity of the Gulf Stream matches well with observations 
and the homogenized pool south of this area loosely 
corresponds with the g minimum described by Keffer (1985). 
This q minimum is also evident in the earlier, hand 
contoured chart of McDowell et al (1982) (Figure 3). 

Figures 8 and 11 display an equatorward decrease in modelled 
potential vorticity vice the poleward decrease in g, north 
of the North Equatorial Current, shown in Figure 3. This 
can be attributed to the limited data set available to 
McDowell et al (1982). The absence of sources of low 
potential vorticity east of this low g region results in 
the water taking on values reflective of the thermocline. 
These values are then subducted and advected southward (note 
the g = 10 contour, Figure 8). The homogenized nature of 
the South Atlantic may be due to the 300 m average depth of 
the density surface, minimizing the influence of surface 
conditions. The homogenization of the south central Indian 
Ocean (Figure 10) can be attributed to this also, with the 
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average depth of the o e = 26.4 surface at 350 m. 

The model Pacific Ocean {Figure 9) and the observed 
Pacific Ocean (Figure 12) both illustrate a dynamic and well 
mixed circulation. Of note in the western North Pacific 
(Figure 9) are the eddy-like features associated with the 
Kuroshio. Their appearance in five-year-averaged data is 
evidence of the strength and permanence of these features in 
the model. The northern and southern tongues of high q 
along the eastern Pacific Ocean are possibly model 
representation of the dynamics of the California Current and 
related upwelling. 

In evidence also are equatorward gradients of q that 
develop on the western side of the subtropical gyres. This 
makes these areas vulnerable to baroclinic instability 
(Keffer, 1983). As with the observations, the gradients 
appear to be due to thermohalrne processes vice the dynamics 
of the recirculating water. 

Notably absent are the isolated lenses of high 
potential vorticity in the eastern North and South Atlantic 
tropical zones. Keffer (1985) attributes these "domes" to 
vertical diffusion (Figure 11). In the model, vertical 
mixing is parameterized according to the method of 
Pacanowski and Philander (1981), a Richardson number 
dependent method for tropical oceans. This parameterization 
may induce smoothing of features smaller than the Richardson 
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The model Indian Ocean (Figure 10) characteristics are 
similar to those of the observed Indian Ocean (Figure 13). 

A feature that may be better resolved by the model is the 
high q region west of Australia. It has been theorized 
that Leeuwin Current water, upon contact with the ACC, 
either mixes with the ACC and continues eastward around 
South Australia or is somehow reflected back into the Indian 
Ocean. Although observed q (Figure 13) does not resolve 
this issue, the flow paths indicated in the model output 
would support the latter. This suggests the need for more 
detailed observations in this area. 

B. UPPER MIDTHERMOCLINE (LAYER B) 

The pressure of the central density surface (a e = 26.75 
in the Atlantic and Indian Oceans, 26.5 in the Pacific 
Ocean) of the upper-midthermocline layer is shown in Figures 
14-16. This layer, theoretically below the deepest 
penetration of wintertime convection, should have no mid¬ 
gyre sources of potential vorticity. Internal processes, 
such as salt fingering may be the exception. Nearly all of 
the recirculation is confined to regions equatorward of the 
outcrops. The recirculating waters are not exposed to 
direct atmospheric forcing and are susceptible to 
homogenization. The surface outcrops of this layer have 
moved poleward. The average maximum depth of this central 
density surface is 600 m. Also evident in each of the 
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basins is the characteristic eastward shoaling of the 
density gradient. Of interest in the eastern North Atlantic 
is the westward bulge of the = 26.75 surface (Figure 14). 
This feature would normally be attributed to Mediterranean 
outflow. The eastern boundary of the model Atlantic Ocean 
is continuous. Therefore, deep thermal and haline restoring 
must account for the missing Mediterranean outflow in the 
model. This feature is not seen in the Keffer (1985) 
analysis of the Levitus (1982) data (Figure 17). This area 
is shown as a maximum in thermal and haline deep restoration 
in the Semtner and Chervin (1992) plots. 

Figures 18-20 illustrate the corresponding charts of 
potential vorticity. All of the oceans show some degree of 
homogenization, with the Pacific Ocean displaying the 
greatest amount. However, the model ocean tends less 
towards homogenization than that observed in the Levitus 
(1982) data (Figures 17 and 21). Again evident in the 
eastern North Atlantic (Figure 18) are the effects of deep 
thermal and haline restoration working to account for the 
missing Mediterranean outflow. 

The tendency toward q homogenization appears to be 
typical. However, the mechanisms for homogenization differ 
between basins. In the North Atlantic, the entire gyre 
shallower than o e = 27.2 is seasonal. That is, winter 
surface densities are generally greater than this value. 
Consequently, the entire subpolar gyre is missing from 
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Keffer's analysis (Figure 17) and the outcrop window lies at 
the front between these two gyres. There is considerable 
downpumping of properties from the outcrop into the gyre 
interior by Ekman pumping. In the model, as in Keffer's 
results, it appears that North Atlantic homogenization is 
the result of an intricate arrangement between the geography 
of the outcrop window and that of the Ekman pumping field. 

By contrast the analogous density surface in the North 
Pacific (Figures 15 and 21) has far less exposure to the 
atmosphere. The outcrop region is confined to a relatively 
small region in the northwest and no simple, direct 
ventilation of this density surface is possible. 

Keffer (1985) attributes the large values of q in the 
subpolar gyre of the North Pacific to a light, fresh layer 
of water in the top 150 m caused by precipitation exceeding 
evaporation. The density layer shown in Figure 21, lies at 
the base of this layer and the f5p/8z values are 
correspondingly high. This precipitation and subsequent 
vertical mixing acts as a strong source of potential 
vorticity to the northern gyre and, consequently, q is not 
homogenized. The corresponding model results. Figure 19, 
show higher values of q in the North Pacific, yet not as 
large as those calculated from observations. The difference 
may be attributable to seasonal and monthly forcing of 
temperature and salinity values in the model. The Levitus 
(1982) fields were spline interpolated to daily values for 
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updating of forcing every three days in the running model. 
This may have introduced additional smoothing and error. 

The high values of q in the northern gyre are not 
efficiently communicated to the subtropical gyre lying to 
the south: fluid particles have their potential vorticity 
reset before they can make an excursion to the adjacent 
gyre. 

Potential vorticity in the South Pacific subtropical 
gyre appears more homogenized in the model (Figure 19) than 
in observations (Figure 21). From these figures it can be 
seen that the area of the winter outcrop of this layer is 
more extensive in the South Pacific than in the North 
Pacific. Figure 19 displays the possible outcome of a 
competition between direct ventilation and a tendency 
towards homogenization in recirculating gyres. Within the 
model South Pacific, the homogenization is complete. In the 
model North Pacific only at the deepest levels, where the 
outcrops no longer surface, is there a tendency towards 
homogenization. 

C. MIDTHERMOCLINE (LAYER C) 

The pressure of the central density surface (<T e = 27.15 
in the Atlantic, Pacific and Indian Oceans) of the 
midthermocline layer is shown in Figures 22-24. This layer 
lies at the base of the wind-driven gyre, in the 
midthermocline (nominal gyre depth of 800 m). It was 
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generally chosen to represent the core of various 
intermediate waters such as the Antarctic Intermediate Water 
or the lighter variety of Subpolar "Mode" Water described by 
McCartney (1982). The anomalous feature in the eastern North 
Atlantic, associated with Mediterranean outflow water, is no 
longer evident at this depth, approximately 700 m. 

Figures 25-27 illustrate the corresponding charts of 
potential vorticity. By this level, the homogenization in 
the North Atlantic (Figure 25) has disappeared and more 
certain signs of direct ventilation reappear again. In the 
observations (Figure 28) the subtropical recirculation of 
Subpolar Mode Water described by McCartney and Talley (1982) 
is evident in the uniform small values of g in the 
northeast corner of the gyre. Once formed in the western 
part of the European basin, presumably south of the zero 
wind-stress curl line, this water is then subducted 
southward between the o e = 27.0 and 27.2 surfaces to appear 
in the gyre as a relative g minimum. Due to the deep 
thermal and haline restoration processes in the model 
eastern North Atlantic, this is not as clearly represented 
(Figure 25). 

The picture in the South Atlantic is strikingly 
different. At first glance, the gyre appears relatively 
well mixed compared to the North Atlantic, the range of g 
values spanning perhaps three-quarters of the total equator- 
to-outcrop range. But a closer examination of observations 
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shows lower g water along 40° S lying poleward of slightly 
higher q water along 30° S - the opposite of the usual |$- 
controlled arrangement. The two waters can be traced back 
to separate source regions. The lower g water appears to 
be McCartney's (1977) “Sub-Antarctic Mode Water" with 
origins in the south-eastern Pacific (see Figure 29), while 
the higher g water appears to originate from the Agulhas 
region (see Figure 28) and can be finally traced back to the 
Indian Ocean (see Figure 30). Keffer's (1985) analysis 
tends to support McCartney's (1977) conclusions on the 
origins of Antarctic Intermediate Water (AAIW). Both 
studies rely on essentially the same technique: tracing low 
g water from the South Pacific to the South Atlantic. 

The picture of this density layer in the model South 
Atlantic is not as clear. The two separate water types may 
be inferred from the shape of the 9.0 g contour (Figure 
25). Accepting this, it does appear that potential 
vorticity properties are being injected or diffused from 
adjacent oceans, rather than a surface outcrop, as in the 
North Atlantic. Note that in the South Atlantic this 
surface outcrops at or slightly south of the zero line of 
wind-stress curl (Baker, 1982), making direct ventilation by 
downward Ekman pumping difficult. However, the location of 
the zero wind-stress curl line is questionable as Southern 
Hemisphere wind field data are relatively unreliable. 

As we have seen, this layer is well ventilated in the 
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North Atlantic and only 'indirectly" ventilated in the South 
Atlantic. However, in the North Pacific (Figure 26) 
isolation is complete and the layer does not outcrop at all. 
Despite its somewhat greater depth than the North Atlantic, 
not only are the North Pacific gyres homogenized, but g is 
homogenized between gyres as well. Both subtropical and 
subpolar gyres equilibrate to a small range of g (see the 
9.0 contour), extending from 25° to 55° N latitude, and 
across the full 10,000 km expanse of the North Pacific. The 
subtropical circulation cannot be driven by the "direct 
ventilation" mechanism. Instead, momentum is transported 
downward via eddy stresses (Rhines and Holland, 1979) where 
it acts to spin up the gyre. Hence, the gyre is indirectly 
driven by the wind. 

Unlike the o e = 26.5-26.7 layer above, at this level 
the South Pacific is nearly homogenized, being filled with 
g = 6 contours except for the intrusion of g = 7 water in 
the southeast. The slight tendency toward nonhomogeneity 
shows itself as g contours that slope from east-southeast 
to northeast. Note the low g values in the southeast 
corner - this is the source of the low g water seen in the 
South Atlantic in Figure 25. 

The Indian Ocean (Figure 27) resembles the South 
Pacific in that the g field is largely homogenized in the 
South with a slight tendency towards northwest to southeast 
ramping (note the g = 9 contour). The small zonal pockets 
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of low q water centered on 43° S (Figure 30) are attributed 
to an artifact of the data base in Keffer's (1985) analysis. 
These are not evident in the model results. 

O. LOWER MIDTHERMOCLINE (LAYER D) 

The pressure of the central density surface (o e = 27.4 
in each of the basins) of the lower-midthermocline layer is 
shown in Figures 31-33. In the North Pacific (Figure 32), 
there is no northern surface outcrop of this layer. The 
density surface itself is relatively flat, with an average 
depth of 1000 meters. The Indian Ocean (Figure 33) has its 
southern outcrop along 50° S, the northern boundary of the 
ACC. The surface's maximum depth is 1400 m in the 
southwest. The surface outcrop in the North Atlantic (Figure 
31) has shrunk and moved northward. The average depth of 
this relatively flat density surface is 900 m. 

Figures 34-36 illustrate the corresponding charts of 
potential vorticity. Figure 34 shows the q field of the 
Atlantic ocean in the o e = 27.3-27.5 layer. This layer was 
chosen because it generally outcrops poleward of the zero 
wind-stress curl line and hence will be exposed to Ekman 
suction (rather than Ekman pumping) within the outcrop zone. 
A still deeper layer might illustrate this better, but may 
add the additional complication of topographic blocking of 
contours by the Mid-Atlantic Ridge. This layer acts as the 
"source" waters for Worthington's (1970) description of the 
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Norwegian Sea as a "mediterranean basin". Warm, light 
surface water flows into the Sea at this density, perhaps 
aided by Ekman suction, and, through atmospheric and 
oceanographic processes is converted into colder, denser 
(and slightly fresher) water to be returned at depth. 

Even at this great depth (over 1000 m) the North 
Pacific g field (Figure 35) is still very homogenous. This 
agrees well with observations (Figure 38). Eddy mixing is 
probably no stronger in the North Pacific than in the North 
Atlantic and may even be weaker due to the much larger 
regions farther from the boundary current. The 
homogenization process is complete due to the high degree of 
recirculation within this layer. 

Unlike the North Pacific, this density layer outcrops 
within the South Atlantic, South Pacific, and Indian Oceans 
(Figures 34-36). Nevertheless, the g field is still very 
nearly homogenous (although the South Atlantic still shows 
signs of influence from adjacent oceans). Despite the 
outcropping, the density layer may be effectively isolated 
from ventilation anyway, either because of the relatively 
weak Ekman suction within the outcrops or because the 
Antarctic Circumpolar Current acts as an effective barrier 
to ventilation by winding g contours around the Antarctic 
continent. The flow path from the eastern North Atlantic to 
deep within the subtropical gyre can be seen in each of the 
Atlantic g charts (Figures 4, 34 and 37). The broad 
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homogenized region in the southern Indian Ocean (Figure 36) 
corresponds well with observations (Figure 39). The east- 
west q contours in the central Atlantic and Pacific Oceans 
and the northern Indian Ocean, seen in Layer C and again in 
Layer D, are characteristic of weakly driven layers that do 
not outcrop. The potential vorticity contours are dominated 
by f = f 0 + Py, indicating a weakly driven flow, with the 
earth's vorticity dominating the potential vorticity fields. 
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VI. CONCLUSIONS 


Charts of potential vorticity (g) within four 
different density layers, calculated from a high resolution 
model that used known bathymetry, observed ship winds, and 
historical temperature and salinity data to simulate global 
ocean circulation were compared with potential vorticity 
calculated from a computerized climatological atlas. 

Overall, the representation of global ocean circulation 
with 0.5° grid size and biharmonic mixing is reasonably 
good. Nevertheless, improved resolution is needed to model 
the instabilities of mid-latitude western boundary currents 
and those of the Antarctic Circumpolar Current. Initial 
results from an improved 1/6° grid size simulation are 
promising. 

New 1/4° and 1/6° simulations, developed at the Naval 
Postgraduate School and the Los Alamos National Laboratory 
in recent months and currently running on the Cray Y-MP/8 at 
the National Center for Atmospheric Research and the 1024- 
node CM-5 at Los Alamos, are showing promising results in 
providing improved resolution of ocean circulation. The 
1/6° model now extends up to 7 8° North and the deep 
diagnostic forcing terms have largely been removed from both 
models. The 1/4° resolution model treats the oceanic 
barotropic mode by the free-surface method of Killworth et 
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al (1991), vice the previous rigid-lid assumption. This 
reduces computer requirements with grid sizes smaller than 
1/2° and allows more realistic geometry with many islands. 
Additionally, the models use observed wind and thermohaline 
forcing together with satellite altimetric data. Analysis 
of model output will show whether these modifications have 
improved representation of potential vorticity fields. 

The density structure of the model is quite similar to 
observations. It might be argued that the model was started 
with the observed density structure and should remain near 
it. However, most of the comparisons made here are in and 
above »"he permanent pycnocline, areas being actively 
modified by surface processes, either directly through 
convection, or indirectly through ventilation and 
subduction. After accounting for the weak restoring below 
710 m and the parameterization of surface heat and moisture 
fluxes, the model results can be considered essentially 
unconstrained. Information in the thermocline range 25-710 
m, outside of polar latitudes, has established itself freely 
over a twenty year simulation (Semtner and Chervin 1992). 

As with the density structure, the potential vorticity 
of the model solution in the density layers above the main 
pycnocline agrees well in the pattern and values with the 
estimates made from observations. The expected 
homogenization occurs below the directly forced layers and a 
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ventilation path from the northeastern Atlantic is evident. 

The deep thermal and haline restoring terms in the 
eastern North Atlantic, there to account for the missing 
Mediterranean outflow, produce an unrealistically large 
signature, when interpreted from the potential density and 
potential vorticity charts. The removal of the deep 
diagnostic forcing terms in a more recent model run may 
solve this problem. Integrating a Mediterranean Sea 
regional model into the global model might also provide more 
realistic results, providing the narrow Strait of Gibraltar 
is treated properly. 

The lack of resolution of the Arctic Ocean causes 
problems which robust diagnostic forcing cannot completely 
remedy. Understanding the effects of the Arctic Ocean and 
its marginal seas on the production of North Atlantic Deep 
Water and the separation of the Gulf Stream, would be 
improved by adding the Arctic domain to the model. As noted 
in Semtner and Chervin (1992), the model Gulf Stream 
separates from the North American continent farther north 
than observed. This may account for the disparity in 
potential vorticity contours and values seen in this area. 
Improved resolution in higher latitudes thus hopefully 
improving treatment of the Labrador current, may improve the 
positioning of the model Gulf Stream. 

Convection and other high latitude processes seem 
poorly treated, even with a seasonal cycle of forcing. More 
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effective treatments of these processes will be required in 
the future. To examine present-day ocean circulation, 
diagnostic forcing can be enqployed to augment convective 
adjustment and maintain deep and bottom waters of the North 
Atlantic and the southern ocean. The results of a 1/4° 
model run with the deep diagnostic forcing terms removed 
showed improvements in resolving deep and bottom water 
currents. 

In the southern high latitudes, however, the lack of 
winter hydrographic data and the unreliable wind stress 
data, especially for a differentiated quantity like wind 
stress curl, make it difficult to evaluate parameterizations 
and model performance. It has come to the point of model 
output exhibiting interesting southern ocean features, not 
previously observed, whose existence are then confirmed with 
further observations. 

Away from boundaries where relative vorticity is known 
to be significant, the charts of potential vorticity appear 
to be good indicators of the dynamics that persist within 
particular density intervals. The transition from the 
direct wind-driven circulation within the gyres to the 
source-sink regime of the main thermocline is observed in 
each of the basins. 

Four separate flow domains are present. In regions of 
Ekman downwelling the surface density distribution is 
prescribed. This regions breaks further into a three-fold 
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subdivision. First, there exists a constantly refreshed, 
ventilated region in which all fluid layers are in motion. 

On the periphery of this ventilated region, shadow zones 
emerge traced by potential vorticity contours emerging from 
meridional boundaries. These shadow zones may bound regions 
of uniform potential vorticity or regions which are 
stagnant. The fourth zone occurs in regions of Ekman 
upwelling. The surface density field there is no longer 
arbitrarily prescribed. Instead, density surfaces are 
pulled, one after another, to the surface in response to 
Ekman suction. 

Although the regions of ventilation and homogenization 
are spatially distinct, they interact and affect one another 
non-linearly since the regions must combine to support the 
Sverdrup transport. 

A valuable feature of using potential vorticity on 
density surfaces to determine flow paths is it allows three 
dimensional flow analysis. With larger vertical density 
variations than in the horizontal throughout the majority of 
the global ocean, vertical resolution is an important 
factor. At 20 layers, the model provides a fair 
representation of vertical thermocline structure. Increased 
vertical resolution may improve that representation 

Further improvements in the physics of the simulation 
and vertical resolution will improve the model's 
representation of oceanic circulation and potential 
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vorticity. In all, the model represents general flow 
characteristics as determined by charts of potential 
vorticity quite well. Recent model upgrades, as mentioned 
above, and planned improvements in forcing and 
parameterizations will continue to improve the model 
representation of ocean circulation. Successive 
improvements come closer to resolving disturbances at the 
local scale of deformation as well as through some portion 
of an inertial cascade, which is one of the main goals of 
ocean modeling. 
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Figure 1. South Atlantic Ocean. Contours of 2ft sin(8a e /5z) 
(DPV/p contours) on the potential density surface o„ = 30.75 gf 1 . 
The potential density surface is referred to a depth of 750 m. ; 
< ^-ff erence of ten units corresponds to a difference of 10' 8 gibin' 1 
s' 1 (From Behringer, 1972). 






Figure 2. South Atlantic Oc 
(DPV/p contours) on the potential 
The potential density surface is 
(From Behringer, 1972). 


ian. Contours of 2£2 sin(8cr e /8z) 
density surface a 9 = 32.10 gf 1 . 
eferred to a depth of 1000 m. 








Figure 3. Potential vorticity (10 -13 cm^s* 1 ) between the a e = 
26.3 and 26.5 surfaces. The bold line represents the outcrop of 
the 26.0 surface during the observation period. The stippled 
area lies between the winter outcrops of the two surfaces. The 
nonlinear control of q contours is clearly evident. Fluid can 
move great distances meridionally, while conserving its potential 
vorticity. Note the northward decrease in q in the North 
Equatorial Current region. This marks this area as a likely 
region of baroclinic instability (From McDowell et a 1, 1982) . 
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Figure 4. Potential vorticity between the <T e = 27.3 and 
27.6 surfaces. Other features as in Figure 3. The q contours 
stretch from northeast to southwest, allowing easy communication 
of fluid between high latitudes and the deep subtropical Atlantic 
(From McDowell efc al, 1982). 
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Figure 11. Potential vorticity of the a 9 = 26.3-26.5 
interval (Layer A) , Atlantic Ocean (10' 13 ends' 1 ) . The striped 
regions mark the winter outcrop window of this layer. Note the 
low g region in the Western North Atlantic. This is the 18° 
Water thermostad of Worthington (1959) (From Keffer, 1985). 
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Figure 12 : Potential vorticity of the o e = 26.05-26.25 interval 
(Layer A), Pacific Ocean. Other features as in Figure 11. Note the high 

q "tongue" entering the eastern sides of the North and South Pacific (From 
Kg t r6r f 1985) • 












Figure 13. Potential vorticity of the (T e = 26.05-26.25 
interval (Layer A), Indian Ocean. Other features as in Figure 11 
(From Keffer, 1985). 
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Figure 17. Potential vorticity of the <J 8 = 26.5-27.0 
interval (Layer B), Atlantic Ocean. Other features as in Figure 
11. Note the regions of homogenized potential vorticity in the 
North and South Atlantic (From Keffer, 1985). 
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Figure 21. Potential vorticity of the o« = 26.25-26.75 interval 
er B), Pacific Ocean. Other features as in Figure 11. Note that while 
subpolar and subtropical gyres of the North Pacific are homogenized, 
e is a sharp contrast of q between them (From Keffer, 1985). 
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Figure 26. Potential vorticity of the o, « 27.0-27.3 interval 

(Layer C), Pacific Ocean. 

















70*N 


100*W 60®W 20*W 2Q°E 



Figure 28. Potential vorticity of the a 9 = 27.0-27.3 
interval {Layer C), Atlantic Ocean. Other features as in Figure 
11. Note the intrusion of low q water into the northeastern 
North Atlantic. This is the recirculation of subpolar mode water 
described by McCartney (1982) (From Keffer, 1985) . 
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Figure 30. Potential vorticity of the = 27.0-27.3 
interval (Layer C), Indian Ocean. Other features as in Figure 11 
(From Keffer, 1985). 
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Figure 35. Potential 
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Potential vorticity of the o # ■ 27.3-27.5 interval 
(Layer D), Indian Ocean. 
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Figure 37. Potential vorticity of the <T e = 27.3-27.5 
interval (Layer D), Atlantic Ocean. The heavy dashed line in the 
North Atlantic is the line of zero wind stress curl. Other 
features as in Figure 11. Note that in the North Atlantic there 
is a tendency for the q contours to run from deep in the 
subtropical gyre into the outcrop region. This facilitates the 
replacement of Norwegian Sea overflow water (From Keffer, 1985). 
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Figure 39. Potential vorticity of the cr 9 = 27.3-27.5 
interval (Layer D) , Indian Ocean. Other features as in Figure 
igoci^°^ e that q is still largely homogenized (From Keffer, 
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